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Abstract. The standard Super-Kamiokande analysis uses an estimator for particle identification by which 
it discriminates e (v e ) from ^(Vfj,). Use of this estimator has led to the claim of a significant deficiency 
of n (u^), suggesting the existence of neutrino oscillations. We investigate three areas of concern for the 
Super-Kamiokande estimator: the separation of the spatial part from the angular part in the probability 
functions, the neglect of fluctuations in the Cherenkov light in different physical processes due to the 
C"*^ ■ charged particles concerned, and the point-like approximation for the emission of Cherenkov light. We 

£\j ' show that the first two factors are important for the consideration of stochastic processes in the generation 

of the Cherenkov light, and that the point-like assumption oversimplifies the estimation of the Cherenkov 
light quantities. We develop a new discrimination procedure for separating electron neutrinos from muon 
neutrinos, based on detailed simulations carried out with GEANT 3.21 and with newly derived mean 
angular distribution functions for the charged particles concerned (muons and electrons/positrons), as 
well as the corresponding functions for the relative fluctuations. These angular distribution functions are 
constructed introducing a "moving point" approximation. The application of our procedure between the 
Oh, discrimination between electron and muon to the analysis of the experimental data in SK will be made in 

a subsequent paper. 

PACS. 13.15.+g Neutrino interactions - 14.60.-z leptons 

x 

■ 1 Introduction GeV. As the results of them, they have found significantly 

different zenith angle distribution between for muon and 
The possible existence of neutrino oscillations is one of the electron, namely, muon deficit and attributed such dis- 
most important issues in particle astrophysics as well as crepancy to the neutrino oscillation between muon and 
elementary particle physics at the present time. Among tau - As the most new one, they give sin 29 > 0.92 and 
the positive and negative results reported for neutrino os- 1-5x10 3 eV < Am 2 < 3.4 x 10 3 eV at 90% confidence 
dilation, experimental results for atmospheric neutrino by level [2]. 

Super Kamiokande (hereafter, we abbreviate simply SK) (3) Also, they have analyzed Upward Through Going 
has special position in the experiments concerned, because Particle Events and Stopping Particle Events. Most phys- 
it is said that they have given the decisive and clear evi- ical events under such category could be regarded as ex- 
dence for the existence of neutrino oscillation. The reasons clusively the muon (neutrino) induced events, not electron 
as follows: (neutrino) induced events, because the effective volume for 

(1) They carried out the calibration experiments for muon is much larger than that for electron due to longer 
the discrimination between muon and electron by electron range of muon irrespective the discrimination procedure 
accelerator beam whose energies are well known and estab- between muon and electron which is indispensable for the 
lished the clear discrimination between muon and electron analysis for Fully Contained Events and Partially Con- 
for SK energy region concerned [T]. tained Events. Also, in this case, they have given the same 

(2) Based on the well established discrimination proce- parameters for neutrino oscillation which are obtained in 
dure between muon and electron, they have analyzed Fully the analysis of Fully Contained Events and Partially Con- 
Contained Events and Partially Contained Events, whose tained Events [2] . 

energies covered from several hundreds MeV to several Through three different kinds of the experiment per- 

formed by SK, all of which are constructed upon the well 

Send offprint requests to: established procedure, it is said that SK has given clear 
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and definite evidence for existence for the neutrino oscil- 
lation. 

The analysis of Fully Contained Events and Partially 
Contained Events is closely and inevitably related to the 
discrimination procedure between electron and muon. Be- 
cause the frequency of muon events with some energy 
occurred inside the detector is nearly the same as that 
of electron events unless ossillation exists and, therefore, 
the precise discrimination procedure between electron and 
muon is absolutely necessary 

Considering the great impact of SK experiment over 
other experiments concerned and theoretical physics, we 
feel we should examine the validities of the experimental 
results performed by SK, because nobody has examined 
them in the most comprehensive way, solely due to char- 
acter of huge experiment, although the partial aspect of 
SK had been examined in fragmental way [3] . 

However, Mitsui et al have examined the validity of 
the discrimination procedure by SK and have pointed out 
the necessity of fluctuation effect into the discrimination 
procedure between muon and electron by SK [4]. 

We have examined validities of all the SK experiment, 
adopting quite different approach from the SK procedure. 

In order to interpret the detected events (i.e. to define 
the CAUSES for the CONSEQUENCES) one should first 
go from CAUSES to CONSEQUENCES, i.e. solve a DI- 
RECT problem, and then back from CONSEQUENCES 
to CAUSES, thus solving an INVERSE problem. SK solve 
the DIRECT problem in a rather simplified way and mostly 
concentrate on the INVERSE problem. We go both ways 
with reasonable care. 

For the purpose, we have performed computer numer- 
ical experiments for examination of the validities of SK 
experimental results, constructing the virtual SK detector 
in the computer, the scale of which is same as the real SK 
detector. 

Concretely speaking, firstly, we have constructed the 
virtual small scale of SK detector, the scale of which is the 
same as the real SK detector (small scale) for being ex- 
posed to the accelerator beam, in the computer, for testing 
the validity of the discrimination between muon and elec- 
trons adopted by SK. SK utlize their results by simple 
exterpolation. 

Based on the imformation obtained from the virtual 
small scale of the detector, we have constructed the full size 
scale of the virtual SK detector, the size of which is the 
same as that of the real SK detector in the computer and 
have really checked whether the clear discrimination which 
SK assert is possible or not. In this case, we take account 
of the difference in the size between the pilot detector for 
the accelerator beam and the real (full size) detector. The 
methodology adopted by us is described in the present 
paper. In a pair of subsequent papers, we really compare 
our results under our methodology (present paper) with 
SK results. Further, the analysis of both Fully Contained 

1 The detailed study for the discrimination between electron 
and muon by using accelerator was carried out in Kasuga [5], [6] 
and Sakai[7]. 



Events and Partially Contained Events and the analysis of 
both Upward Through Going Muon Events and Stopping 
Muon Events in the numerical computer experiments will 
appear in the subsequent papers. 

An analysis of Fully Contained Events and Partially 
Contained Events is given in a subsequent paper. An anal- 
ysis of Upward Through Going Muon Events and Stopping 
Muon Event will be presented elsewhere. 

From the methodological point of view for the analysis 
of the SK data, there should be two essential differences 
between the SK procedure and our procedure in the gen- 
eral analysis of physical events. 

The first: as we already mentioned above, SK oversim- 
plifies the solution of the DIRECT problem, namely, their 
mean models of both electron and muon events are too far 
from reality to be used in the INVERSE problem solution, 
while our mean models are reasonably accurate. The sec- 
ond: SK neglects fluctuation effects in physical processes 
in which the Cherenkov light is produced, while we con- 
sider them as correctly as possible, which we discuss in 
the section 2 and the subsequent sections. 

In this paper, we limit our discussion to the discrim- 
ination procedure between muon and electron in the SK 
experiment (the item(l) in the Introduction) 

The contens of the paper are organized as follows: In 
the section 2, we examine the SK standard discrimination 
procedure in detail. We examine the point-like approxi- 
mation in the electron shower due to electron neutrino 
adopted by SK and conclude that such approximation 
leads to serious error in the energy determination of the 
electron neutrino. Next, we examine the probability func- 
tions for the pattern part and the angular part in both 
muon and electron and conclude that such separation of 
the pattern part from the angular part is not adequate 
and further, there are, a priori, no reasons why the prob- 
ability functions for muon obey the same type of that for 
the electron, even if such separation is valid. 

There are many factors which produce various errors 
in the SK discrimination procedure. The simple agebraic 
sum of such errors not always denote right resultant errors. 
In our discrimination procedure, we examine the resultant 
errors from the SK standard discrimination procedure, for 
example, the error in the vertex points of the neutrino 
interactions and their directions and more precise values. 

In the section 3, we develope more suitable discrimina- 
tion procedure instead of the SK standard procedure and, 
for the purpose, costruct the mean angular distribution 
functions for Cherenkov light due to muon and electron 
and the corresponding relative fluctuation functions which 
make it possible to estimate the degree of the separation 
between muon and electron, errors of the vertex points for 
muon-like event and electron like events, the error of the 
directions due to these events. In a subsequent paper, we 
develop the general procedure for estimating the concrete 
errors for the neutrino events concerned and give finally 
various errors for the neutrino events quantitatively by 
comparing our resluts with SK results. 
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2 Examination of the SK discrimination 
procedure between muon and electron 



The approach for the discrimination between muon and 
electron adopted by SK is as follows. 

(1) SK calculate Cherenkov images for electron and muon 
events by statistical and deterministic numerical methods 
but while constructing the pattern recognition procedure 
all fluctuation data and some important features of the 
mean models are completely ignored which are discussed 
below. Such manipulation leads to serious error in both 
pattern recognition between electron and muon and their 
energy estimation. 

(2) Based on the Cherenkov light of the particle con- 
cerned mentioned above, SK constructed the probabil- 
ity function, ESTIMATOR for the particle identification, 
which is composed of pattern part and the angular part 
in the separated form. 

The probability function for pattern, P pa ttern(e(iJ,)) , 
which will appear in the section 2.2, denotes the pattern of 
the events concerned, namely, spatial and angular image of 
the events concerned. However, the variety of the images 
of the events concerned in SK discrimination procedure 
exclusively come from the fluctuations in the photoelec- 
trons in the PMT, but not from the physical processes of 
the events concerned which produce the Cherenkov light, 
while in our procedure the image of the events concerned 
are essentially governed by the physical proccesses con- 
cerned which are completely neglected in the SK proce- 
dure. 

The probability function for the angular part is in- 
troduced into the total probability (estimator for particle 
identification), being separated from that for pattern part. 
This is unnatural, however, because the pattern is of the 
concept on the space-angular structure. And it is natu- 
ral that the concept of the angular should be included in 
the concept of the pattern. In our procedure, we construct 
the concept of the pattern as the angular-space image in 
which the spatial one is interrelated with the angular one 
(See, section 3.3). 

Further, SK adopt the same type of the probability 
functions for muon and electron. However, there are no 
reasons why the probability function for muon obeys that 
for electron, because the generation mechanism of the to- 
tal Cherenkov light is quite different in muon from in elec- 
tron. 

In present paper, we would emphasize the neglect of the 
fluctuation effects and the oversimplification in the prob- 
lems adopted by SK lead to serious errors on the discrim- 
ination between muon and electron, which are shown in 
the present section and subsequent sections. 



2.1 The Cherenkov light calculation for electron and 
muon. 

2.1.1 Examination on the Treatment of Cherenkov light for 
electron by SK 

The SK analysis approximates the sources of the Cherenkov 
light due to electron showers as being point. In the calcu- 
lation of the source for the Cherenkov light due to shower 
particles, the SK analyze electron events in average val- 
ues, neglecting fluctuation and further assigns them to 
be point-like source. Such a treatment is not valid for the 
analysis of the real physical events concerned in the SK de- 
tector, because the sources for the Cherenkov light due to 
shower particles definitely have geometrical extents. The 
geometrical extent of electron cascade shower as a light 
source is not negligible compared with the scale of the SK 
detector. 

Thus, the introduction of the point-like approxima- 
tion into the construction of the ESTIMATOR for parti- 
cle identification generally invites not negligible error into 
the analysis of experimental events. In what follows, we 
examine the limitations of the point-like approximation 
by SK. Here, we restrict our examination to the direct 
Cherenkov light, because the scattered Cherenkov light is 
not essential for the examination of the logical structure 
of the SK discrimination between electron and muon and 
we are exclusively interested in the logical structure of the 
SK procedure. 

In order to construct the ESTIMATOR for electron 
identification, the SK calculate the mean number of pho- 
toelectrons which are produced by the Cherenkov light 
due to an electron-primary cascade shower (electron neu- 
trino) using the Monte Carlo method. Thus, SK obtains, 
Ni^ exp (direct), the expected number of photoelcctron from 
the direct Cherenkov light initiated by an electron in the 
i th PMT at a distance and an angle, in the following : 

N it e Xp (direct) = 
= a e ■ N M C(e hPe ) ■ ■ exp • /(0),(1) 

where NMc(6i,Pe) denotes the mean number of the pho- 
toelectrons received in a circular area of 50 cm diameter 
located on a sphere of 16.9 meter in radius by the Full 
Monte Carlo simulation, a e is the normalization factor, li 
is the distance from the particle position to the i th PMT, 
0i is the angle of the i th PMT from the particle direc- 
tion, and f{0) is the effective photo-sensitive area of the 
PMT. Here, L is the attenuation length of the Cherenkov 
light in the SK detector, which is taken 100 meter. The 
third term in the right-hand side of Eq.(l) shows the at- 
tenuation of light and its power index is estimated from 
Monte Carl simulation. The value of 7 is 2.0 in Takita 
(p83, in [9]) and Kasuga(p32 in [5]) and Sakai(p41 in [7]), 
and 1.5 in Kasuga(p71 in [B]). A numerical value of 2.0 is 
used in older analysis while 1.6 is adopted in more recent 
work. The expression of Eq.(l) is an oversimplification on 
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Fig. 1. The dependence of the Cherenkov photon density due 
to the electron shower on the distance for differnt angle from 
the direction of an incident 1 GeV electron with the point-like 
approximation, normalized with the values at 16.9m in Fig- 
ure 2. The angles considerd are 5, 42 and 85 degree. Sampling 
numbers of Monte Carlo simulations are 10000 per each angle. 
See the text for detail. 
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Fig. 2. The dependence of the Cherenkov photon density due 
to the electron shower on the depth for differnt angle from the 
direction of an incident 1 GeV electron without the point-like 
approximation. The angles considerd are 5, 42 and 85 degree. 
Sampling numbers of Monte Carlo simulations are 10000 per 
each angle. See the text for detail. 



the estimation of photoelectrons by the Cherenkov light 
when we consider the real behaviour of the electron cas- 
cade shower. In order to clarify the oversimplification of 
the SK procedure, let us compare it with our procedure 
in the level of Cherenkov light, but not in the photoelec- 
tron level, because it is enough for the purpose concerned. 
Here, we coulud define the following expression which cor- 
respond to Eq.(l), 

N exp ,cher,app(dir€Ct) 
= a' e ■ N M C,cher{0,Pe) X X ex P ' ^ 

where NMC,cher(6,Pe) is the mean Cherenkov photon den- 
sity which correspond to NMc(@i,Pe) in Eq.(l), a' e is an- 
other normalization factor and other parameters are same 
as in Eq.(l). Namely, Eq.(2) gives the Cherenkov photon 
density for the electron cascade shower under the point- 
like approximation which is utlized by SK. The depen- 
dence of the Cherenkov photon density N exp ^ c her.app{direct) 
on the distance I is essentially determined by (16. 9/l)~< 
rather than exp(-^/i), because the attenuation of the Che- 
renkov light is small (100 meter of attenuation length in 
SK). Thus, the Cherenkov light depends strongly on the 
inverse power of the distance in the SK procedure. 

In Figure 1, the dependence of the Cherenkov photon 
density due to the electron cascade shower on the distance 
(/) are given for differnt angle from the direction of an 
incident 1 GeV electron with the point-like approximation 
using Eq.(2). The energy value of 1 GeV is a typical energy 
in the SK detector. 

In order to examine the validity of the point-like ap- 
proximation by SK, we have carried out the calculation 
of N eX p, c her,FMC {direct), the Cherenkov photon density 
without the point-like approximation in Eq.(3) for the 



comparison with Eq. (2) . Instead of the point-like approx- 
imation, we simulate the electron cascade shower exactly 
in the stochastic manner, which we call [without point- 
like approximation]. Then, the Cherenkov photon density 
generated by the shower elctron is calculated by help of 
GEANT3.21. 

Namely, Eq.(3) denote the straightforward expression 
for the Cherenkov light density which is directly derived 
from the electron cascade shower without the approxima- 
tion of (16.9/Z) 7 and exp(-//i). Thus, the definition of [ 
without point-like approximation] denotes that we simu- 
late electron cascade shower and accompanied Cherenkov 
light density in three-dimensional way at any distance 
and angle as exactly as possible. Namely, we calculate the 
whole structure of the Cherenkov photon density, taking 
into account the attenuation of light, as well as that of the 
electron cascade shower. Then, we obtain 

cher,FMc(direct) — NMC\cher,FMc(0,Pe,l), (3) 

where NMC,cher,FMC'(6,Pe, I) denotes the Cherenkov pho- 
ton density at the distance / without introducing the point- 
like approximation, which is obtained by the Monte Carlo 
method, taking into account the attenuation length of the 
Cherenkov light correctly. 

In Figure 2, the dependence of the Cherenkov photon 
density on the distance without point-like approximation 
is given for the different angles from the direction of the 
incident elctron using Eq.(3). In [without point-like ap- 
proximation], the dependence on the distance and angle 
and the attenuation are automatically included into func- 
tions themselves. 

From the Figure 1, it is clear that the dependences of 
the Cherenkov photon density on the distance in the SK 
are same irrespective of the angles, which simply reflects 
the separation of the distance part from the angular part 
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Fig. 3. Cherenkov photon density as a function of angle for Fig. 4. Cherenkov photon density as a function of angle with- 

different distances with the point-like approximation normal- out the point-like approximation for 1 GeV primary electron, 

ized with the values at 16.9 m in Figure 4 for 1 GeV primary Other parametersare the same as in Figure 3. 
electron. 



in Eq.(2). Comparing the Figure 2 with Figure 1, it is eas- 
ily understood that in the region beyond about 10 meter, 
the point like approximation is well approximation as far 
as the tendency in concerned, except the absolute values, 
while in the region inside about 10 meter, this approxima- 
tion does not hold anymore, due to the actual longitudinal 
and lateral structure of the electron cascade showers. In 
particular, for direction within 5 degree, the effect of the 
lateral spread of the electron cascade shower as well as 
that of the definite longitudinal length becomes effective 
and the actual dependence on the distance is deviated 
largely from that in the point-like approximation. 

In Figure 3 and Figure 4, we give the angular de- 
pendence of the Cherenkov photon density for different 
distances [with point like approximation] (Eq. (2)) and [ 
without point like approximation ](Eq.(3)), respectively. 
In Figure 3, the shapes of their angular dependence is ex- 
actly same irrespective of their distance, because their an- 
gular part is separated from their distance part, as shown 
in Eq.(2). Such separation of the angular part from the dis- 
tance part does not reflect real situation of the Cherenkov 
light density for electron, as shown in Figure 4. 

In Figure 4, we give the corresponding ones to Figure 3 
in the case of without point-like approximation. It should 
be noticed that for the region of smaller than 5 m the ten- 
dency of the angular dependence on angle is largely devi- 
ated from that for the region of larger distances. This re- 
flects the fact that the Cherenkov photon density is largely 
influenced by the longitudinal and lateral structure of the 
real electron cascade shower. From figure 1 to 4 normal- 
ization are made at 16.9 meter. 

In Table 1, we give the ratio of the Cherenkov pho- 
ton density with the point like approximation to that 
without the point-like approximation for the exponent of 
1.5(Kasuga,p74[6i]). In Table 2, we give the same quanti- 
ties for the exponent 2.0, the older value used in the SK 
analysis (Takita,p83. 9], Sakai,p41[7j). The ratios are nor- 



Table 1. The ratio of the Cherenkov photon density from 
electron showers with the under point-like approximaton 
(Eq.(2),7 = 1.5) to corresponding ones without the point- 
like approximation (Eq.(3)). The primary enegy of electron is 
1 GeV. 



angle 






distasnce 




(m) 


(degree) 


5.00 


10.00 


16.90 


20.00 


40.00 


5 


0.175 


0.553 


1.000 


0.948 


1.452 


10 


0.191 


0.548 


1.000 


0.949 


1.457 


20 


0.200 


0.525 


1.000 


0.956 


1.494 


30 


0.231 


0.495 


1.000 


0.973 


1.583 


40 


0.549 


0.640 


1.000 


0.939 


1.505 


50 


0.757 


0.756 


1.000 


0.887 


1.135 


60 


0.715 


0.739 


1.000 


0.893 


1.166 


70 


0.708 


0.736 


1.000 


0.893 


1.171 


80 


0.708 


0.737 


1.000 


0.892 


1.167 


85 


0.706 


0.737 


1.000 


0.891 


1.164 



Table 2. The ratio of the Cherenkov photon density from elec- 
tron showers with the point-like approximaton (Eq.(2),7 = 2.0) 
to corresponding ones without the point-like approximation 
(Eq.(3)). The primary enegy of electron is 1 GeV. 



angle 






distasnce 




(m) 


(degree) 


5.00 


10.00 


16.90 


20.00 


40.00 


5 


0.322 


0.719 


1.000 


0.872 


0.944 


10 


0.352 


0.712 


1.000 


0.872 


0.947 


20 


0.368 


0.682 


1.000 


0.879 


0.971 


30 


0.425 


0.643 


1.000 


0.894 


1.029 


40 


1.009 


0.832 


1.000 


0.863 


0.978 


50 


1.393 


0.983 


1.000 


0.815 


0.738 


60 


1.314 


0.960 


1.000 


0.821 


0.758 


70 


1.301 


0.956 


1.000 


0.821 


0.761 


80 


1.301 


0.958 


1.000 


0.820 


0.758 


85 


1.299 


0.958 


1.000 


0.819 


0.757 
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malized to the values in Eq.(3) at the 16.9 meter in the 
both cases. 

Comparing Table 1 and Table 2, it is clearly under- 
stood that the ratios depend strongly on the values of 7. 
If the point like approximation is valid, then, the ratios 
in Table 1 and Table 2 should remain around 1.0 with 
satisfactory allowance for every distance and every an- 
gle. However, it is clear from the tables such situation is 
never realized. SK guarantee the accuracy of the energy 
determination for the event is within ±2.6 % (Ishitsuka 
±2 % (p29 in [TD]), Kameda ±2.5 % (p62 in [IT]), Oku- 
mura ±2.5 % (p45 in [12]), Messier ±2.5 % (p92 in [13]) 
and Kasuga ±2.6 % (p53 in [6])). This is quite far from 
the reality. In conclusion, Eq.(l), the expression for the 
photoelectron utilized by SK, does not consider the longi- 
tudinal and lateral structure of the electron cascade and, 
therefore, an oversimplification which could not estimate 
the photoelectrons by the Cherenkov light correctly. 



2.1.2 Cherenkov light for muon 

In contrast to the case of electron, SK calculates the Che- 
renkov light, taking into account the definite extent of 
muon range, namely, without point-like approximation. 
However, still, SK neglect the fluctuation effect. Here, let 
us examine the validity of the expression adopted by the 
standard SK analysis for a muon event, in which the ex- 
pected number of photoelectrons in the i th PMT pro- 
duced by a muon is expressed as Takita(p84 in [S]), 

Sakai(p41 in [7]), Kasuga(p74 in [6]), Kibayashi(p71 in 

i). 



N i>exp (direct) 



k(sin6i + k x (£)) 



x sin 2 6i + Ni,knock(0i 



xexp 



L 



x /(©), 



(4) 



where a M is the normalization factor. The second term 
in the right-hand side originates from the ionization en- 
ergy loss dE/dx in water. Taking into account the change 
in the photon density which is caused by the change in, 
corresponding to the energy loss, dxsinO + ld9, expresses 
the intensity variation of the Cherenkov photons. 

Ni,knock(9i) shows the number of expected photons 
from knock-on electrons as a function of, which is esti- 
mated by a Monte Carlo simulation. It is easily understood 
from Eq.(4) that the logic for producing photoelectron in 
the muon is definitely deterministic. In other words, the 
SK analysis neglects fluctuation effects on the generation 
of the Cherenkov light from muon completely. However, 
the muon losses its energy and changes the direction in 
stochastic way, namely, knock-on, the decay-product elec- 
tron and multiple scattering. Therefore, we could not ne- 
glect fluctuations for the muon event. As one example of 
muon behavior, we show the range fluctuations for muon 
in Figure 5. The neglect of fluctuations for muon causes 
uncertainty in the estimation of the muon energy. 



300 * 5GeV/c \i 



60 



20 



10 




1600 1700 1800 1900 2000 2100 2200 2300 2400 

H range, cm 

Fig. 5. Range fluctuation of 5 GeV muons. The number of 
simulation is 300. 



2.2 Pattern recognition procedure adopted by SK 

In previous subsection, we examine accuracies on the Che- 
renkov light quantities for electron and muon, which is 
related to the estimation of the energies of the particles 
concerned directly for the Fully Contained Events. In the 
present section, we should like to say these quantities are 
of major importance for the type definition and geometry 
recognition, both in SK and our case. Here, we examine 
pattern recognition procedure for the particle concernrd. 



2.2.1 Principle of the ESTIMATOR of particle identification 
adpted by SK 

SK defines the probability, Probi(N exp , N b s ), suitable for 
the particle identification. It is the probability to observe 
the number of photoelectron, N b s , for an expected mean 
number of photoelectron, N exp . The N t s actually devi- 
ates from the ideal Poisson distribution due to fluctuations 
in amplification processes of the PMTs. The probability 
function is simply given by: 



Prob 



1 



2ttct 



exp 



(N obs - N exp )' 
2a 2 



(5) 



where Probi is the probability function for the i th PMT. 
Thus, SK define the likelihood functions for e-like event 
and mu-like event is the following: 



L e = J] Fro5 *(e), 
0i<(i.5xe c ) 



(6) 



9«<(1.5x0 c ) 

where Prob^e) is calculated assuming the event is due 
to electron. And Probi(p) due to muon. In order to com- 
bine the information of the Cherenkov opening angle with 
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the information of the ring pattern, L e and L M are trans- 
formed into x 2 distribution so that SK obtain the follow- 
ing functions for spatial part and angular part. 



2/„\ ,,2/,,M \ 2 



1 ( X 0) - min[ X 2 (e),x J (fj)] 



Ppattern(e) = exp < - 



Ppattern(p) = exp • 



Pangie(e) = constant x exp < —- I — 



(7) 



Pangie (/■*) = constant x exp • 



1 f 8expW — ^obs 
A0 



(8) 



Combining (6) with (7), SK obtains probability functions 
for electron and muon finally in the following: 



P(e) = Ppatternie) X P a „gle{e), 
Pit 1 ) = Ppatternin) X Pangle(^) ■ 



(9) 



A ring is more mu-like than e-like, if > P{e) and 
vice versa. 



2.2.2 Examination of the components of the ESTIMATOR 
adopted by SK 

Let us examine the validity of Eq.(5) in the SK proce- 
dure. [1] They neglect fluctuation effects in physical pro- 
cesses for electron and muon. SK introduce Eq.(5) for the 
detection of the particle concerned. However, they con- 
sider only fluctuation effect coming from the amplifica- 
tion processes in the PMT, namely, electro-mechanical ef- 
fect of the hardware detector, and never any fluctuation 
effect coming from the physical process for electron and 
muon, namely, fluctuation from shower particles for elec- 
tron (neutrino) and range straggling of muon for muon 
(neutrino) which play essential role for formation of vari- 
ety of the pattern onto the detector for electron and muon. 
Namely, SK describe the behavior of electron and muon in 
the mean values. Consequently, varieties of the Cherenkov 
light pattern due to electron and muon are produced only 
through "Poison distribution in the amplification process 
in the PMT, the position of the neutrino interaction and 
the direction of the incident neutrino. However, essential 
character of the varieties of the Cherenkov light due to 



2 In our opinion, a priori, there are no reasons why does 
the probability function for muon obey the same type of the 
probability function for electron. 



electron and muon are only produced through the fluctu- 
ation in the physical processes of the particle concerned. 
If wc consider fluctuation effect correctly, Eq.(l) defined 
by SK should be replaced by 



P e (N : Eo) = 
= Pe(N e , obs: N e , exp (E )) x P e (N,N ei0bs ), (10) 

where P e (N : Eq) denotes the probability for an electron 
with energy Eq to produce a cascade shower in which JV 
photoelectron are produced by the Cherenkov light due 
to shower particles. P e (N ey0bs , N e ^ exp {Eo)) is the proba- 
bility for an electron with primary energy Eq and the 
mean number of the Cherenkov photon, N e-exp (Eo), to 
produce the observed number of the Cherenkov photons, 
N ei0 b s and P e {N, N ey0bs ) denotes the probability for N^ obs , 
number of the Cherenkov photons to produce finally the 
number of the photo-electron N. P e (N e ^ obs , N ei< , xp (E )) 
shows the fluctuation effect in a electron casade shower 
and P e (N,N e ^ obs ) correspond to Eq.(5) in SK procedure. 
For muons, the expression which correspond to Eq.(8) is 
given as follows: 



P^(N : E ) = 
= Pp,(N ll+et0hs ,N ll+e , exp {Eo)) x P p {N,N p+e , obs ){ll) 

where Pp,(N : Eq) is the probability for a muon of pri- 
mary energy Eq to produce the number of the photo- 
electron N by the Cherenkov light due to the muon and 
its accompanying electron (knock on electron and others) . 
Pp,(N p+e ^ obs , N p+e ^ exp (E a )) shows the fluctuation effect in 
the physical processes of muon, namely, the range strug- 
gling of muon. Corresponding to Eq.(5), P e (N, N ey0 b s ) and 
Pp(N, Np+e^obs) may be approximated by Gausian distri- 
bution with different N. However, P e {N e ^ obs ,N e ^ exp (Eo)) 
should be different from P A1 (A r A1+e!0bs , N M+eiexp (Eo)), be- 
cause the physical process for producing the Cherenkov 
light for electron is clearly different from that for muon. 
Consequently, there are no reasons, a priori, why 
Pe{N e ,obs, N eiexp (E )) and P p (N p+e , obs , N p+etexp (E )) 
should obey the same Gaussian distribution. In conclu- 
sion, Eq.(9) lacks in theoretical background as the proba- 
bility function in the practical application to the analysis 
of experimental data. 



2.2.3 Resultant errors introduced by SK discrimination 
procedure 

In previous subsection, we examine errors introduced into 
the Cherenkov light quantities for electron and muon. The 
final purpose of the construction of the ESTIMATOR 
adopted by SK is [a] the discrimination between electron 
and muon, [b] decision of the vertex position in the neu- 
trino event and decision of the direction of the individual 
neutrino event and [c] the decision of the momentum of 
the particle concerned. Therefore, one may be very inter- 
ested in the resultant errors which SK procedure invite 



8 



V.I. Galkin et al.: A Theory of Discrimination in SK 



due to inadequency of the ESTIMATOR, the neglect of 
the fluctuation effect and assumption of the point-like ap- 
proximation. 

The answer is as follows: we could not estimate indi- 
vidual error in the SK procedure, for example, the error 
due to the the oversimplification in the treatment of the 
electron cascade shower. The individual error is intermin- 
gled with each other so that we get the resultant errors 
which are never the simple sum of the individual error. 
Namely, we could get the resultant errors only, for exam- 
ple, the error of the energy measurement in the neutrino 
events, the errors on the vertex in the neurino event, the 
error on the direction of the neutrino even and so on. 

In the next section, we develop our rigourous procedure 
which does not include inconsistency of the logic. 

3 Principle of our discrimination procedure 
between electron and muon: Application of 
pattern recognition theory to e-mu 
discrimination. 

3.1 Overcome to the defect of the SK procedure 

In the previous section, we clarify that the estimation of 
the Chrenkov light quantities due to the electron event 
and the muon event as well as the discrimination of their 
pattern are not reliable. Because, both electron events and 
muon events are recognized without taking their fluctua- 
tion effect which may kill the real variety of their pattern. 
In the electron events, the point-like approximation may 
lead serious misestimation of their as shown in Table 1 
and Table 2. In order to overcome the defects inherent 
in the SK discrimination procesure, we develop our dis- 
crimination procedure as exaxtly as possible, taking the 
stochastic characters in the physical processes concerned 
carefully into account. 

Here, we restrict the development of our theory at the 
level of the Cherenkov light, but not at the photoelectron 
level which are directly related to the actual experimental 
result as in the SK procedure. The reasons are as follows: 
The most important of our present paper is to clarify the 
essential defect adopted by SK and propose the direction 
of the possible alternative procedure by which SK proce- 
dure could be replaced. 

In this case, it is enough for us to examine physical 
quantities in the level of the Cherenkov light leaving out 
the consideration of photon conversion into photoelectron. 
Such simplification reduces the overall uncertainties of pa- 
rameter estimates and, thus, our results give the lower lim- 
its for the errors in the analysis of the neutrino events con- 
cerned. To illustrate the scale of uncertainty added by the 
light-to-electron conversion, in Figure 6 we give the proba- 
bility distributions of the number of the Cherenkov photon 
for given number of photoelectron. Of course, while apply- 
ing our methods to the analysis of the real experiment one 
should take photoelectron fluctuations into account. 
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Fig. 6. fig:6 Probability distributions of the number of the 
Cherenkov photons for given numbers of photoelectron. The 
number of simulation is 100. 

3.2 Construction of a suitable procedure for the 
discrimination of muons from electrons 

The examination described in the previous section leads us 
to conclude that the probability functions and the approx- 
imations used in the standard SK analysis are not reliable 
for the discrimination of muons from electrons. Therefore, 
a more reliable procedure, constructed on a solid theoreti- 
cal basement, is required. Here, we will outline a more re- 
liable and suitable algorithm for discrimination of muons 
from the electrons. First, we simulate the physical pro- 
cesses for the Cherenkov light as exactly as possible using 
GEANT 3.21, and construct new mean angular structure 
functions for the Cherenkov light for pattern recognition 
of the Cherenkov light due to both the muon and electron. 
Second, we derive their relative fluctuation functions for 
the Cherenkov light distribution which are critical for con- 
structing an estimator for particle identification. 

In our mean angular structure functions, we calculate 
the Cherenkov light due to both muon and electrons along 
the direction of the the primary muon and electrons. Then, 
followings arc considered. 
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(1) As for the mean longitudinal size of /i-track/clectron 
shower, we adopt the total longitudinal length (T e defined 
lately) in which 99.5 percent of the total Cherenkov light 
are considered. 

(2) The total longitudinal length along which the particles 
concerend emit the Cherenkov light is emitted is devided 
into a number of eual segments. The number of the di- 
vision is dependent on the particle type, primary energy 
and the required accuracy of the image presentation. 

(3) We calculate the mean angular distribution, F^'^(9), 
and its relative fluctuation S^'^(0) for each segment i,and 



following to construct the angular distribution function 
for the total Cherenkov light. Here, we construct the mean 
angular distribution function for the Cherenkov light for 
electron. 

In a large water tank detector which is large enough 
for the dimensions of the electron showers concerned, we 
calculate the development of total Cherenkov light due to 
the electron shower. 

We exactly simulate individual angular distributions 
for the Cherenkov light using a combination of GEANT 3.21 
with the calculation tools developed by us. To construct 



construct the mean angular distribution functions, F e ^{6, t) the mean angular distribution function for the Cherenkov 



and thcri relative fluctuations, <5 e,Al (6*, t), as function of ra- 
diation angle, 9, and water layer thickness, t. The item(3) 
will be discussed in a subsequent paper 

We neglect the lateral distributions of particles in both 
e and \i events, i.e., we assume that all Cherenkov light 
is emitted at the event axis. This makes our approach in- 
valid for events close and approximately parallel to the 
walls of the detector, but it is adequate for the particle 
discrimination quality estimates. The much simpler ap- 
proach adopted in the standard SK analysis would result 
in larger errors in type and geometry reconstruction pro- 
cedures, particularly in case of peripheral events when the 
longitudinal development of e-shower//z-track Cherenkov 
light angular distribution is most prominent. 



light and its relative fluctuation function accurately, we 
need to simulate a large number of showers, 10000 to 
20000. Let T e be the total length of an electron shower 
initiated by a primary electron. We divide T e in segments, 
the length of which depends on the primary energy of the 
charged particle concerned and which is taken as 40 cm in 
present calculation. 

Shower particles are produced according to the ex- 
act simulation procedure of the electron cascade shower, 
but in our models of the Cherenkov light, both the mean 
aogular distribution function and corresponding angular 
distribution relative fluctuations function are attributed 
to a certain segment. When calculating the pattern im- 
ages with the help of the models, all Cherenkov photons 
generating within the segment are thought to start from 
Note that angular distribution functions for the CherenkoYhe middle of the segment. This approximation is called 



light are universal functions for a given particle type and 
energy and can be used for calculating mean pattern im- 
ages and their deviations for any required event geometry 
in any water tank. 

Further, it should be noticed that we need not introduce 
Ppattem(e(fJ-)) and P a ngie(e(fJ-)) in the procedure for the 
discrimination as in Eqs.(7) and (8) in the final probability 
function, and directly get the final probability function. 



3.3 The construction of the mean angular distribution 
functions of Cherenkov light for muons and electrons 



the "moving-point" (or "multi-point") approximation, in 
contrast to the single point approximation adopted in the 
standard SK analysis. 

We simulate the cascade shower initiated by the pri- 
mary electron and calculate the Cherenkov light from shower 
particles in each segment and obtain the angular distribu- 
tion function for the Cherenkov light by distributing all 
the Cherenkov photons emitted from segment k over the 
bins of a histogram N e (9, E , k) in 9 to estimate: 
The Mean angular distribution function 



From the view point of pattern recognition, we need more 
accurate angular distribution functions for the total Cherenkov 
light due to both muons and electrons, rigorously taking 
into account fluctuations inherent in both types of event, 
to confirm the discrimination between muons and elec- 
trons. For this purpose, we construct angular structure 
functions for the Cherenkov light in the following way. 



F e {9i,E ,k) 



(N e (9 l ,E ,k)) 



(12) 



where (...) denotes the average over a large event sam- 
ple, 9i is the center of mass of the i-th histogram bin, and 
AQi is the solid angle of the z-th bin. 
The Relative fluctuation function for the angular 
distribution 



3.3.1 The construction of the mean angular distribution 
function due to a primary electron and its relative 
fluctuation function. 

An electron due to electron neutrino interaction produces 
an electron shower in which shower electrons are distributed 
over some range in space, each of which produces Cherenkov 
light. It is not suitable for such a range to be approximated 
by the point-like as assumed in SK, on which we have al- 
ready clarified. Instead of the SK point-like approxima- 
tion, we introduce a 'moving-point' approximation in the 



5 e (9 t ,E ,k) = 



yJ(N2{0i,Eo,k)) - (N e (6i,E ,k)y 



(N e (9 t ,E ,k)} 



(13) 

Then we fit the mean angular distribution (12) for each 
segment separately by a function: 



F e (9,E ,k) = 

= iq{[A + B3/(i+B-\8-C\( b4+1 



i-\e-c\) 



)]/[l + B2.0]} 



(14) 
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0.5 GeV e, segment No.1 



fit: F(6)=10**((A + B3/(1 + B*abs(6-C) 
"(B4 + B1*abs(6-C)))) / (1 + B2*6)) 

A=2,2794 
B=0.56371 
C=4 1.6869 
B1 =0.001 275 
B2=0.002186 
B3=4.2706 
B4=0.58105 
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Fig. 7. fig: 7a The angular structure function for 500 MeV 
electrons from segment No.l. See the text for details. 



where the numerical values of A, B, C, Bl, B2, B3 and 
B4 depend on the primary energy of the electron and the 
corresponding segment. 

In Figures 7 to 11, we give the mean angular distri- 
bution function for the total Cherenkov light from each 
segment for a 500 MeV electron. We give the mean an- 
gular distribution function for the total Cherenkov light 
emitted in segments No.l (the interval from to 40 cm 
from the starting point of cascade shower), No. 2 (from 
40 cm to 80 cm from the starting point), No. 3 (80 cm to 
120 cm), No.4 (120 cm to 160 cm), and No.5 (160 cm to 
200 cm), in Figures 7, 8, 9, 10 and 11, respectively. 

The electron shower produced by a 500 MeV electron 
is rather small and the essential character of the cas- 
cade shower is determined in segments 1 to 5. In Fig- 
ure 7 (segment No.l), which corresponds to the initial 
stage of shower development, shower particles are ener- 
getic and near the core of shower axis and so give rise to a 
strong peak around the Cherenkov angle due to energetic 
particles initiating the cascade shower. In Figure 8 (seg- 
ment No. 2), the cascade shower reaches shower maximum, 
where number of shower particles reaches a maximum, 
and the total Cherenkov light produced in this segment is 
roughly the same as that in segment No.l. However, the 
average energies of shower particles in this segment are 
smaller than those in segment No.l, and they are more 
scattered due to multiple scattering and so there is more 
deviation from the original Cherenkov angle. 

The situation in segment No. 3 (Figure 9) is roughly 
the same as that in segment No.2(Fig.8). It should be 
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0.5 GeV e, segment No.2 

o - fit: F(6)=10**((A+ B3/(1 + B*abs(6-C) 
**(B4 + B1*abs(8-C)))) / (1 + B2*6)) 

A=2.7373 
B=0.0371 
C=41.9232 
B1=1.22e-4 
B2=0.00207 
B3=1 .972 
B4=1.1718 
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Fig. 8. fig:7b The angular structure function for 500 MeV 
electrons from segment No.2. See the text for details. 
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0.5 GeV e, segment No.3 



o - fit: F(6)=10"((A + B3/(1 + B*abs(6-C) 
"(B4 + B1*abs(6-C)))) / (1 + B2*e)) 

A=2.416 
B=0.0271 
C=42.057 
B1=3.06e-4 
B2=0.00217 
B3=1 .920 
B4=1 .0981 




10 20 30 40 50 



60 70 80 

9, deg 



90 



Fig. 9. fig:7c The angular structure function for 500 MeV 
electrons from segment No.3. See the text for details. 
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0.5 GeV e, segment No.4 

1n 5[ o - fit: F(6)=10**((A + B3/(1 + B*abs(6-C) 
1 u ; "(B4 + B1 *abs(6-C)))) / (1 + B2*6)) 



C=42.050 
B1=4.91e-4 




10 r 
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e, deg 

Fig. 10. fig:7d The angular structure function for 500 MeV 
electrons from segment No.4. See the text for details. 

0.5 GeV e, segment No.5 

s[ o - fit: F(e)=10**((A + B3/(1 + B*abs(9-C) 
1 u : **(B4 + B1 *abs(6-C)))) / (1 + B2*e)) 

= A= 1.986 

B=0.0254 
C=42.127 
B1=2.96e-4 




10 r 
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Fig. 11. fig:7e The angular structure function for 500 MeV 
electrons from segment No.5. See the text for details. 



emphasized that the amounts of Cherenkov light emit- 
ted in segments No. 2 and No. 3 are greater than that in 
the segment No.l (near the starting point of the cas- 
cade shower) at all angular regions except the original 
Cherenkov angle. Even in the segments No.4 (Figure 10) 
and No.5(Figure 11), where electron shower is attenuated 
rapidly, the contributions of Cherenkov light from outside 
the original Cherenkov angle could not be neglected com- 
pared with that from the segment No.l, which is too near 
the starting point of the cascade shower. In conclusion, 
as the Cherenkov light produced in every segment con- 
tributes to the edge of the Cherenkov ring, we cannot say 
that the Cherenkov light at the edge of the Cherenkov ring 
comes exclusively from the starting point of the cascade 
shower (the vertex point of the electron neutrino reac- 
tion), which is the assumption adopted in the SK analysis 
(See Figure 2.9 in Sakai,[7]). Here, we have paid particu- 
lar attention to a 500 MeV electron, however, such char- 
acteristics hold irrespective of the primary energy of the 
electron. 

In order to calculate the relative fluctuation function 
for the angular distribution accurately, a large sample 
of simulated events is needed: we used 20000 events for 
E <500 MeV and 10000 events for E >500 MeV. 

In Figures 12 to 16 we construct the relative fluctu- 
ation function for the Cherenkov light from a 500 MeV 
electron. The relative fluctuation functions are extremely 
important for the pattern recognition of electrons and 
muons. Therefore, we should determine them as accurate 
as possible, which was why we simulated such a large num- 
ber of events. 

In Figures 12, 13, 14, 15 and 16 we give the relative 
fluctuations for the Cherenkov light from segments 1, 2, 
3, 4 and 5, respectively. 

In Figure 12 (segment No.l), the shower particles are 
relatively energetic so that the relative fluctuation is a 
minimum near 42°, the original Cherenkov angle, as ex- 
pected. In Figure 13 (segment No. 2), the situation is the 
same as in Figure 12. It is interesting to note that the 
minimum disappears in the segment No.3(Fig. 14) and 
a maximum appears at the original Cherenkov angle in 
segments No.4(Fig.l5) and No.5(Fig.l6). However, glob- 
ally speaking, the angular dependence of the fluctuation 
is weak and it is relatively similar at all angles. Of course, 
this does not mean that shower particles are produced 
isotropically, losing the direction of the initiating particle. 

The fluctuations in the electron events which are de- 
scribed above are quite different from those in muon events 
(see next subsection). 

We need to take these relative fluctuations into ac- 
count for more reliable discrimination between electrons 
and muons. This is the reason why we calculate the rel- 
ative fluctuations after simulating each energy 10000 to 
20000 times. 
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Fig. 12. fig:8a The relative fluctuations in the angular dis- Fig. 14. fig:8c The relative fluctuations in the angular dis- 
tribution of the Cherenkov light for 500 MeV electrons from tribution of the Cherenkov light for 500 MeV electrons from 
segment No.l. See the text for details. segment No. 3. See the text for details. 
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Fig. 13. fig:8b The relative fluctuations in the angular dis- Fig. 15. fig:8d The relative fluctuations in the angular dis- 
tribution of the Cherenkov light for 500 MeV electrons from tribution of the Cherenkov light for 500 MeV electrons from 
segment No. 2. See the text for details. segment No. 4. See the text for details. 
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Fig. 16. fig:8e The relative fluctuations in the angular dis- 
tribution of the Cherenkov light for 500 MeV electrons from 
segment No. 5. See the text for details. 
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o - fit: F(e)=10**(A*exp(-B*(9-C) 2 )) 
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Fig. 17. fig:9a The angular structure function for 500 MeV 
muons from segment No.l. See the text for details. 



3.3.2 The construction of the mean angular distribution 
function due to a primary muon and its relative fluctuation 
function. 

In the same way as in the preceding sub-section, we obtain 
the mean angular distribution function F^{9, Eq, k) for the 
Cherenkov light emitted due to the primary muon and its 
relative fluctuation 5^(9, E ,k). 

We then fit the mean angular distribution for each seg- 
ment in turn by a function: 

F^e,E ,k) = iQ{Ae X p[-B(e-cf]} + 10 [m/(i+B2^)] ^ 

(15) 

where the numerical values of A, B, C, Bl and B2 are 
dependent on both Eq, the primary energy of muon, and 
fc, the segment. 

In Figures 17 to 21 we construct the mean angular 
distribution function for the total Cherenkov light from 
each segment for a 500 MeV muon. We give the mean an- 
gular distribution function for the total Cherenkov light 
due to segments No.l (the interval from to 40 cm from 
the starting point of cascade shower) (Fig. 17), No. 2 (40 cm 
to 80 cm) (Fig. 18), No.3 (80 cm to 120cm) (Fig. 19), No.4 
(120 cm to 160 cm) (Fig.20) and No.5 (160 cm to 200 cm) 
(Fig. 21) Comparing Figure 17 with Figure 7, one can eas- 
ily understand that the peak around 42° in the case of the 
muon is very sharp compared with that in the case of the 
electron. 

This means that the larger contribution of the Cherenkov 
light is due to muon itself and the contributions to either 
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Fig. 18. fig:9b The angular structure function for 500 MeV 
muons from segment No.2. See the text for details. 
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Fig. 19. fig:9c The angular structure function for 500 MeV 
muons from segment No.3. See the text for details. 
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Fig. 21. fig:9e The angular structure function for 500 MeV 
muons from segment No.5. See the text for details. 



0.5 GeV ji, segment No.4 

ln 5[ o - fit:F(e)=10"(A*exp(-B*(9-C) 2 )) 
1U : +10**(B1/(1+B2*6 4 )) 
= A=3.792 

B=0.00220 
C=32.634 
B1 =2.649 
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Fig. 20. fig:9d The angular structure function for 500 MeV 
muons from segment No.4. See the text for details. 



side of the peak are from knock-on electrons and others 
with lower energies in the case of a muon. In the case of the 
electron initiated shower, the cascade of electrons with dif- 
ferent energies and directions produces a broader spread of 
Cherenkov light. As the muon proceeds it loses energy, so 
that the angular distribution for the Cherenkov light be- 
comes wider by multiple scattering while its mean value 
decreases due to the reduction of the Cherenkov angle. We 
can see this in Figures 18, 19 and 20, where the angular 
distributions change greatly between segments, while the 
production of Cherenkov light is maintained, which is a 
different characteristic from the corresponding curves in 
the case of the electron. In Figure 21, we find the disap- 
pearance of the peak of the angular distribution which was 
evident in the preceding segments. The uniformity of the 
distribution in Figure 21 shows that the decay electrons 
are distributed isotropically in this segment due to the 
nature of the three-body decay. According to the calcula- 
tions on transition curves for the differential and integral 
Chcrenkon photons, it is easily understood that there is 
no muon in segment No.5(Fig.21) (160 cm to 200 cm from 
the starting point of the muon). 

Also, we could not infer the starting point of the muon 
(the vertex point of muon neutrino interaction) from the 
edge of the Cherenkov ring, as adopted by the SK group, 
just as for the electron event. Comparing Figure 17 with 
Figures 18,19,20 the contributions from the segments No. 2 
(Fig. 18), No.3(Fig.l9) and No.4(Fig.20) near the edge of 
the Cherenkov ring are comparable with that from the 
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Fig. 22. fig:10a The relative fluctuations in the angular distri- 
bution of the Cherenkov light for 500 MeV muons from segment 
No.l. See the text for details. 



segment No.l. In conclusion, we can determine the vertex 
point of the neutrino interaction, or the direction of the 
incident neutrino, based on the total pattern of Cherenkov 
light only, and not solely from the edge of the Cherenkov 
ring. 

We can see that there is a clear difference in the relative 
fluctuations between the muon and electrons as well as 
that in the angular distribution functions. In Figures 22 
to 25, there are peaks in the relative fluctuation either side 
of the local minimum at 42°. More exactly speaking, the 
strong concentration of the Cherenkov light around 42°, 
which is the consequence of the large contribution by the 
muon gives a deep minimum to fluctuation. The two peaks 
appear as the result of this, and the local minimum shifts 
to smaller angles, accompanied by the two peaks, as the 
muon proceeds, while there is only one peak in the case of 
an electron-initiated cascade. 

However, in Figure 26 where the muon decays com- 
pletely in the segment No. 5 (160 cm to 200 cm from the 
starting point of the muon), the (decay product) electron 
is, on average, distributed uniformly and the distribution 
of the Cherenkov becomes uniform, so that the relative 
fluctuation also becomes uniform. 

Comparing Figure 22 to 26 with Figure 12 to 16, it is 
easily understood that we could expect larger fluctuation 
in muon events than in electron events. The reasons are as 
follows: In electron events, the shower particles produced 
compensate for the effects of fluctuations, while in the 
muon events the single muon concerned bears the effect 
of the fluctuation exclusively due to multiple scattering, 
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Fig. 23. fig:10b The relative fluctuations in the angular distri- 
bution of the Cherenkov light for 500 MeV muons from segment 
No. 2. See the text for details. 
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Fig. 24. fig: 10c The relative fluctuations in the angular distri- 
bution of the Cherenkov light for 500 MeV muons from segment 
No. 3. See the text for details. 
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Fig. 25. fig:10d The relative fluctuations in the angular distri- 
bution of the Cherenkov light for 500 MeV muons from segment 
No. 4. See the text for details. 
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Fig. 26. fig:10e The relative fluctuations in the angular distri- 
bution of the Cherenkov light for 500 MeV muons from segment 
No. 5. See the text for details. 



which are not smeared out. As a result, we could expect 
a larger uncertainty in muon events for both the vertex 
position and the direction of the neutrino than in electron 
events. We discuss this problem in a subsequent paper. 
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